Spin and Current Correlation Functions in the d-density Wave State of the Cuprates 
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We calculate the spin-spin and current-current correlation functions in states exhibiting d x 2_ y 2- 
density wave (DDW) order, d x 2_ y 2 superconducting order (DSC), or both types of order. The 
spin-spin correlation functions in a state with both DDW and DSC order and in a state with 
DDW order alone, respectively, illuminate the resonant peak seen in the superconducting state 
of the underdoped cuprates and the corresponding feature seen in the pseudogap regime. The 
current-current correlation function in a state with both DDW and DSC order evinces a superfluid 
density with doping dependence which is consistent with that of the underdoped cuprates. These 
calculations strengthen the identification of the pseudogap with DDW order and of the underdoped 
cuprates with a state with both DDW and DSC order. 



PACS numbers: 



I. INTRODUCTION 



Recently, a new order termed d-density wave (DDW) 
was proposed for the underdoped regime of the high-T c 
cuprate superconductors §]-||. It was argued that the 
observed anomalies in the pseudogap gap phase are due 
to the DDW and its competition with the rf-wave super- 
conductor (DSC). Recent elastic neutron scattering j|, 
/iSR ||, and NMR/NQR || experiments appear to have 
found some direct evidence for the existence of DDW or- 
der in YBCO by observing the magnetic fields associated 
with this order parameter - in other words, by measuring 
the order parameter itself. In this paper, we wish to elu- 
cidate the nature of the competition between the DDW 
and DSC order paramters by explicitly computing the 
spin-spin and the current-current correlation functions 
in two space-dimensions in the Hartree-Fock approxima- 
tion. Our aim is to provide complementary, indirect evi- 
dence in favor of the DDW scenario by showing in some 
detail that it provides a consistent explanation of the 
existence and behavior of the resonant peak seen in in- 
elastic neutron scattering and also of the evolution of the 
superfluid density with doping and temperature, at low 
temperatures. 

Since we deal with broken symmetry states with order 
parameters, it is possible to argue that their descriptions 
arc adequately described by the Hartree-Fock approxi- 
mation, at least deep within the ordered phases at zero 
temperature. The Hartree-Fock approximation by itself 
cannot be adequate if collective modes are important, 
but these modes could be described in terms of small 
oscillations about the ordered state by suitable random 
phase approximations (RPA). However, critical fluctua- 
tions close to a quantum critical point cannot be treated 
by either Hartree-Fock or RPA. Nonetheless, it is use- 
ful to explore the consequences of the simplest possible 
Hartree-Fock and RPA theories to see if any robust fea- 
tures arise due to the existence of the new order param- 
eter, DDW, and its competition with DSC. 



For underdoped bilayer superconductors Q, the dy- 
namic spin-spin correlation function is peaked above T c 
at the in-plane wavevector Q = (n, 7r) (in this paper we 
shall take the lattice constant to be unity and will also 
set h = c = ks = 1) in both even and odd channels 
with respect to the layers in the bilayer complex. Be- 
low T c , there is a resonant peak at the same wavevec- 
tor and approximately the same energy in the odd chan- 
nel. A resonant peak is also found below T c in optimal 
and overdoped bilayer high-T c superconductors. There 
are many approximate theoretical calculations of the dy- 
namic structure factor Here we concentrate only on 
the aspects germane to our discussion of DDW. We wish 
to stress that this analysis is simply a consistency check 
of the DDW; we do not claim that our explanation of the 
resonant peak in inelastic neutron scattering is unique or 
better than the others. The real test of the DDW hy- 
pothesis is the direct observation of the order parameter 
referred to earlier We find that the experimental 

observations of the neutron dynamic structure factor are 
consistent with the hypothesis that the neutron scatter- 
ing peak - and, hence, the pseudogap - is due to DDW 
order, while the resonance peak in the underdoped super- 
conducting state is due to the coexistence of DDW and 
DSC orders j|. The suppression of the resonant peak by 
a magnetic field perpendicular to the plane suggests that 
superconducting pairing is important for its formation 

The behavior of the superfluid density in the under- 
doped regime of the high-T c materials is another inter- 
esting physical quantity. Many experiments |l0| have 
indicated a rapid collapse of the zero-temperature super- 
fluid density as the doping is decreased below optimal. 
The zero-temperature penetration depth, for example, 
grows rapidly in the pseudogap regime, correlating with 
the suppression of T c with underdoping, yet saturates 
at overdoping in a way reminiscent of a traditional BCS 
superconductor Jll| . This, as we will show, can be ex- 
plained by the competition between DDW and DSC or- 
ders. At finite temperatures the superfluid density is sup- 
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pressed from its zero temperature value at that doping 
fraction. This suppression is linear in temperature (T) 
at asymptotically low temperatures, with a slope that 
is independent of doping. This behavior is captured in 
our calculations by the thermal excitation of nodal quasi- 
particles into the upper quasiparticle band. For heavily 
underdoped samples, there is an intermediate low tem- 
perature regime in which we predict this suppression to 
be proportional to \T. There is some experimental evi- 
dence for this behavior jf2|m YBa2Cu40s, although the 
explanation of these authors involves a proximity model 
of alternating stacked superconducting and normal lay- 
ers. Our analysis is simpler and follows from the nodal 
exciations in the mixed DDW and DSC state. This pre- 
diction can be tested in future experiments in the regime 
in which pseudogap dominates the DSC gap. 

II. ORDER PARAMETER 

Our Hartree-Fock analysis of the physical quantities 
merely requires us to specify a mean field Hamiltonian 
with the proposed order parameters. The actual micro- 
scopic Hartree-Fock analyses to obtain these order pa- 
rameters have been discussed in the past on numerous 
occasions and need not be repeated here; for a recent set 
of calculations, see Rcfs. [||, and references therein. 
Hamiltonians with short range repulsion and superex- 
change have the DDW ordered state (both the singlet and 
the triplet variety) as one of their many possible saddle 
points. If one includes correlated hopping terms, then 
within a reasonable range of parameters, the DDW sad- 
dle point can be stabilized over the other possibilities. In 
the present paper, we shall parametrize the order param- 
eters phenomenologically, and our main conclusions are 
independent of any microscopic self-consistent Hartree- 
Fock equations. 

The singlet DDW state is defined by the order param- 
eter @ 

( c L c k+ q/ 3> = i-^(cosk x - cos k y )5 a p, (1) 

where Cka is the fermion destruction operator for 
wavevector k and spin-index a, and q is the ordering 
wavevector. This order parameter is a generalization of 
the familiar charge density wave (CDW) order parameter 
to the case of angular momentum two, in exact analogy 
to the generalization of the BCS s-wave superconducting 
order parameter to its <i-wave (DSC) counterpart. 

For q = Q, the case of interest to us, the underly- 
ing bipartite square lattice band structure is equivalent 
under the transformation Q > — Q and this forces the 
DDW order parameter to be imaginary [pj. Thus, time- 
reversal symmetry is broken (i.e. the system exhibits 
magnetism) , and the ground state has an array of bond 



currents, which alternate in direction (clockwise, anti- 
clockwise) in the neighboring plaquettes of the 2D bi- 
partite square lattice. The corresponding s-wave CDW 
order parameter is, of course, real. 

III. SPIN DYNAMICS 
A. Pure DDW state 

The Hartree-Fock DDW Hamiltonian, using (1), is 

k(7 kcr 

(2) 

where the DDW gap is given by 
W 

Hk = — (cosfer - cos k y ), (3) 

£k = —2t(cosk x + cos k v ) gives the band structure, and 
ji is the chemical potential. A more realistic band struc- 
ture would include the effect of next-neighbor hopping 
t' . This does not affect our results here, so we drop it for 
simplicity. At half filling, the chemical potential /i = 0, 
while fi takes non-zero negative values as we introduce 
holes into the system. 

At fi = the zero-temperature spin-spin correlation 
function S, at momentum transfer q = Q, is given by 

S(Q, u,n = 0) = 3tt £ S L - 2^Jw£ + eA , (4) 

kerbz ^ ' 

where the integration is over the magnetic or reduced 
Brillouin zone (rbz). This function is peaked at an energy 
equal to twice the maximum value of the gap, 2Wq. As 
we dope holes into the system, at nonzero values of /i, 
Eq. H is changed to 

kerbz ^ ' 

xeL+^/wi + eA (5) 

The only effect of fj, is the depletion of the integration 
region. Hence, as long as \/i\ < Wo, the peak at q = Q 
exists and stays at the same energy 2Wq. 

Note that S(uj) also shows a peak for momentum trans- 
fer q = (n, 0) [Q at the Hartree-Fock level. From a re- 
lation analogous to that of Eq. g, we can see that the 
peak-energy for q = (ir, 0) scales with t. This is much 
higher in energy than the (ir, ir) peak; it is likely to have 
an enormous width (since there is an enormous amount 
of phase space into which it can decay) and be unob- 
servable. This is important because no peak has been 
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observed at q = (tt, 0) in the spin-fluctuation spectrum 
of the pseudogap regime of the high-T c materials. Early 
mean-field decouplings of the Heisenberg model jl4j led 
to an effective t ~ J, so that the peaks at q = (tt, 7t) and 
q = (n, 0) were comparable in energy, in contradiction to 
the experiments. 



B. Pure DSC state 

Below the superconducting transition temperature T c , 
a DSC state is defined by the gap 



Al 



Ac 



cos k x — cos k y ) 



(6) 



corresponding to the DSC order parameter. This order 
parameter, as stressed in reference (jj, can compete and 
coexist with the singlet DDW order parameter. In fact, 
at half filling, the system can be rotated continuously 
from a pure DSC order to a pure singlet DDW order and 
vice versa, without ever having to close the quasi-particle 
gap . (The symmetry between the s-wave counterparts 
of these two types of orders is exact in the negative U 
Hubbard model @.) 

Using (6), the mean- field one-body DSC Hamiltonian 

is 

#DSC = XI ( £k ~ ^) C L C kff + A k c kT C Ul + h ' C - ( 7 ) 
ka k 

The spin-spin correlation function, at fi = and momen- 
tum transfer q, is 



S(q,u,n = 0) = ^ ^ 



1 - 



e k£k4 



kErbz 

X S(UJ - E k - Sk+q), 



+q 



(8) 



where E^ = yf + e£. For q = Q the coherence factor 
equals two and Eq. || becomes identical to Eq. ||. Hence, 
just like the DDW, at \i — the DSC spin-spin correlator 
also shows a peak at twice the gap maximum, 2Ao. 

But for a DSC, at non-zero values of n, S(q, is 
given by 



<5(q,w,/i) 



3tt 



E 



i 



£k£k4 



A k A k , 



kerbz 

X S(U -E k - Sk+q), 



£ki?k+q 



(9) 



where E^ now incorporates \i in the usual way by ek — > 
?k = £k — M- Equation ^ shows a peak at q = Q that 
is shifted in energy from 2Ao. The peak-energy shift is 
small for small , but can be easily checked to scale 
with 2\\i\ as becomes large. Hence, one can clearly 
see that the peak-energy behaves differently with fj, for 
DDW and DSC. 



These Hartree-Fock calculations for the resonance peak 
in the normal state becom e interesting only if the chem- 
ical potential remains pinned close to fi = in the un- 
derdoped regime. (This is not true for the calculation of 
the superfluid density reported in the following section, 
which is insensitive to such assumptions for pi.) There 
is some evidence that this is the case in the photoemis- 
sion measurements of Ino et al. |T^] in La2- K Sr :E Cu04. 
They find that while the chemical potential shift is large 
in the overdoped samples, it is largely suppressed in the 
underdoped regime. Numerical studies of the 2D Hub- 
bard model flj],[l8| also suggest that the chemical poten- 
tial does not shift much for small doping- fractions. In 
the Monte Carlo study performed in reference [0, the 
calculated shift of [i follows A/i = /i cx — x 2 for small 
values of doping, x. The data of reference [|l6| suggest 
that the chemical potential varies with x in the over- 
doped regime as one would expect for a Fermi liquid, but 
not in the underdoped regime. Ideas involving charge 
ordered stripe states [|l9) are suggestive of the segrega- 
tion of doped holes in the boundaries of antiferromag- 
netic domains, thus pinning the chemical potential. The 
presence of charge inhomogeneity, such as stripes [ p0[ , 
or impurity-disorder may pin the chemical potential in 
other multilayer high-T c cuprate systems as well. 



C. Coexisting DDW and DSC 




FIG. 1. The high-Tc cuprate phase diagram proposed in 
Ref. 1. Tjv denotes the 3D- antiferromagnetic transition tem- 
perature. T c and T* are the DSC and pseudogap, or DDW, 
transition temperatures, respectively. In the text (Sec. V), 
the behavior of the spin-s pin correlator is discussed along the 
three paths shown here by the dashed lines. 

From the order parameter competition picture pro- 
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posed in reference M (See, Fig. |l|), one notes that to the 
right of x ~ 0.05 and below T c , the two orders, singlet 
DDW and DSC, coexist and compete for the same re- 
gions on the Fermi surface. The mean-field Hamiltonian 
for the system in the folded Nambu basis plf is 



(10) 



kerbz 



where = (cj^, c^q^, c_ k j_, c_ k -Qj.) and the matrix 
A k is 



peak is located at + A§ for fx = 0. For small fx 

the peak-energy is shifted to higher values, initially by 
a small amount, but finally scaling as 2\fi\ for large val- 
ues of \fx\. We should mention here that one can in fact 
control the behavior of the peak-energy with fx by adjust- 
ing the relative strengths of the order parameters. When 
the singlet DDW order overshadows the DSC order, the 
peak-energy will tend to be pinned to the total order 
parameter magnitude 2^Wq + A§, a behavior charac- 
teristic of the singlet DDW. 



/ (e k - fx) iW k 



A k \ 

-iW k -(e k + /i) -A k 

A k -Ok-M) iW k 

\ -A k ~iW k (e k + fx) J 



(11) 

The matrix A k has four eigenvalues, ±E±(k.) and 
±£ 2 (k), where £i(k) = [( + W£ - fx) 2 + Aft 1 / 3 and 



A k ]V2. 



The spin-spin correlation function is obtained by first 
evaluating the imaginary-time ordered correlator (u> m is 
the bosonic Matsubara frequency) 



(T T S(q,', 



S(-q, -iu) m )) = 



— J2 Tr [G(k,z^)G(k + q,^„ + ^ m )] (12) 

n,k£rbz 



D. RPA Corrections 

If we add a reduced Coulomb repulsion U, which is 
assumed to be renormalized due to particle-particle cor- 
relations, we can drive the system toward antiferromag- 
netism. To model this effect, we have to go beyond the 
Hartree-Fock, and we use a crude RPA form jl| for the 
susceptibility x(q, w) given by, 



x(q»w) 



Xo(q,w) 



i - Uxo{q,u) 



(14) 



to describe the spin dynamics in the system. Here 
Xo(q, oS) is the Hartree-Fock susceptibility. Extracting 
<S(q, lo, fx) from x(q, u, fx), in Fig. || we plot the results as 
a function of frequency for five different values of U. 



where S denotes the spin operator and the G's are the 
4x4 matrix Green functions computed from the Hamil- 
tonian in Eq. [Io|. After analytically continuing the fre- 
quency to the real axis, we extract the imaginary part 
and use the fluctuation-dissipation theorem to get the 
spin-spin correlation function. 

The zero-temperature spin-spin correlation function 
thus obtained, at non-zero fx, for momentum transfer 
q = Q, is given by 



S(Q,uj,fx) = y 

k6rbz 



A k 



x 5 {oj ■ 



£i(k)£7 2 (k) 
E^k) - E 2 (k)) 



(13) 




When Wfc = 0, the coherence factor in Eq. clearly 
matches with the coherence factor in Eq. ^, where we 
have to put q = Q, and the two expressions become 
identical. On the other hand, for A k = 0, one notes that 
£i(k) 



■M 



and£ 2 (k) - 



Keeping in mind that \x is negative in our hole-doped 
system, we note that for \A k + W£ + fx > the coher- 
ence factor in Eq. O is two, while the other choice for 



v4 + 



■ fx < reduces the coherence factor to zero. 



Thus the expression in Eq. |13| becomes identical to the 
expression in Eq. [| for A k = 0. 

One can easily check that S(Q,ui, fx), as a function of 
w, has a peak in the spin-spin correlation function . The 



FIG. 2. iS(Q,OJ,/i) as a function of lo after U is introduced 
into the system with coexisting singlet DDW and DSC orders 
withi = 0.5 eV, A - 02 eV, W =.02 eV and fx set to -.01 eV. 
As U is increased, the structures are peaked at lower and lower 
energies; from right to left U_ = 0.2, 0.217, 0.233, 0.256, 0.284. 
Finally they evolve, as U = U c = 0.32eV is approached, into 
a divergence at zero energy. 

The first thing one notices is that after the introduc- 
tion of U, the correlation function S, which was logarith- 
mically singular when derived from Xo(q, alone, now 
becomes broad. Also, as U is progressively increased, 
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the spectral-weight appears at smaller and smaller en- 
ergies. For U large enough, the correlator starts peak- 
ing up again and finally evolves into a divergence at 
zero energy as U c is approached. All of this behavior is 
qualitatively consistent with what happens to the spin- 
fluctuation spectrum as the doping is reduced in the un- 
derdoped regime of the cuprates. As the hole density 
is reduced, the repulsive Coulomb energy becomes more 
and more important and gradually approaches the crit- 
ical value U c for which the peak in the spin-spin corre- 
lation function evolves into a divergence at zero energy, 
signalling antiferromagnetic order in the system. 

This analysis applies both below and above T c in the 
high-T c cuprate materials. Above the superconducting 
transition temperature T c , Ak will go to zero, but Wk will 
survive up to the pseudogap temperature T*. Hence, in 
the spin-spin correlator the DDW will continue to main- 
tain its signature intact in the form of a peak at twice 
the maximum of its gap. 



The paramagnetic current response is obtained from 
the relation 

A xx (q,iu m ) = {T T j x (q,iuj m )j x (-q,-iu; m )) 
At 



— sin k x sin(k + q) x 

n,k£rbz 

Tr[G(k, iu n )MG(k + q, iuj n + iu 

A W2 

sin k x sin(k + q) x 



f3 ^ 



n,k£rbz 

Tr[G(k, iu n )NG(k + q, iu) n + ico m )N] (18) 

In Eq. [jlj, the G's are the same matrix Green functions 
used in the calculation of the spin-spin correlation func- 
tion and M and N are the matrices 



M = 



/ 1 \ 
0-100 
10 

\ -1 J 



IV. SUPERFLUID DENSITY 

The anomalous behavior of the superfluid density in 
the underdoped regime of the high-T c cuprates has been 
discussed extensively. Here, we address it in the con- 
text of the competition of the order parameters that 
we have been discussing p2|. We will make a phe- 
nomenological assumption about the doping dependence 
of the DDW and single particle gaps, respectively, Wk 
and \J A 2 . + W£ . Following P3[ , we hypothesize that 
the DDW order parameter vanishes above a critical dop- 
ing x c (about 0.2) and follows the formula Wojks — 
aJ*(l—x/x c ) where J* = 980K, x is the doping-fraction, 
and a ~ 0.6 is a constant arbitrarily chosen to set a scale. 
In other words, 



W ~ .049(1 - x/x c ) 



(15) 



in eV; in Fig. we test another model in which Wo cx 

a-x/x c ) 2 . 

We also posit that the single-particle gap remains con- 
stant at its value at x = .05, where superconductivity 
starts developing; then 



W 2 + Ag = $2 = (.037) 2 



(16) 



We calculate the superfluid stiffness for such a model 
using the formula [E4j 

— ^ = (~k x ) - A xx (q x = 0, q y — > 0, iuj m = 0) (17) 
m 

which gives the response to a transverse vector poten- 
tial A x (t). Here, n s is the superfluid density and m* is 
the effective mass, (k x ) is the kinetic energy per site per 
lattice dimension, and K xx {(\^iuj m ) is the paramagnetic 
current-current correlation function. 



N = 



I 10 \ 
-10 
-1 

\ o oi o / 



The superfluid stiffness (and therefore the inverse 
square penetration depth) at T = 0, evaluated in the 
entire doping range, is given by 



m* 



= -32t 2 sin 2 h x W£ ( 1 



kerbz 
1 



el + Al + WZ 



/r 



£a(k)£ 2 (k) 



[£a(k) + £: 2 (k)] : 



-(1 + eW/O) 



S^lim^j^oo sin 2 k x 

k£rbz 

/#(#i(k)) , df(E 2 (k)) 



V rf-Bi(k) 
4t cosfc^ 

k£rbz 

E 1 (k)+E 2 (k) 



dE 2 (k) 



(l + 0(W 2 /t 2 )) 



A 2 . 



W 2 -f 



£i(k)£ 2 (k) 
(l + 0(W 2 /t 2 )) 



(19) 



where E\ (k) and E 2 (k) are the two energy- values given 
before. The first two terms in Eq. [l9] come from the 
paramagnetic current response, and the last term is the 
kinetic energy (or the diamagnetic term). At finite, but 
low temperatures, one can show that 



—AT)-—* 

m m 



(0) « 8t 2 Y sin2 k x 



kerbz 



W 2 
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The 0(WQ/t 2 ) terms derive from the modification of 
both the kinetic energy and current operators in the low- 
energy effective (Hartree-Fock) Hamiltonian associated 
with DDW order. The modified current operator gives 
rise to the matrix N in ( |i~8| ) . Since these terms are small, 
0(Wq /t 2 ), we neglect them in the following discussion. 

Equation [l9] has several interesting properties. First, 
for Ak = Ao ^ and Wk = Wo, that is, for compet- 
ing CDW and s-wave superconductivity, the terms in- 
volving the Fermi functions do not contribute at zero- 
temperature. At /i — 0, Ei(k) = i?2(k), and the same 
coherence factor in the other two terms becomes two. 
Then, by a partial integration one can show that (de- 
noting the first term by K para and the third term by 



j^'dia'j ^para _ _ 



K 



total 



^-dia 



and the full kernel is given by 

K dmM _ W*s when the twQ 

orders are of s-wave type [H , K tot&x yields a superfluid 
density which is maximum when Wo=0 (i.e, Ao = $0) 
and zero when W =$ (i.e, A o =0) p5[ |. 

Returning to d-wave order parameters, at Ao = 0, 
Eq. |l9] reduces to the superfluid density of the DDW 
state, which is, of course, zero. The second term, which 
contains the derivatives of the Fermi functions is crucial 
for the cancellation in this case. For Ao and /1 finite, 
the second term does not contribute at zero temperature 
(except for a contribution from a single point in k-space; 
but the finite and almost constant contribution from this 
is ignored here for simplicity). The system now acquires 
a finite superfluid density. 



The superfluid density derived from Eq. 19 can be cast 
as a function of a single variable, x, by expressing the 
three parameters Wo, Ao and /j, in terms of the doping 
fraction. For the dependence of [i on x, we use [i ~ — x 
(in units in which 2t = 1) as in reference Ji~7[ , but other 
reasonable dependences should give similar results (c/. 
the discussion in the previous section). Hence, the re- 
sults of this section are valid and relevant for the cuprates 
even if \x is not pinned close to zero in the underdoped 
regime. For Wo and Ao we choose the relations given 
in Eq. [l5| and Eq. [16] only because they have received 
some recent experimental support p3[ . But it's worth- 
emphasizing that the qualitative features that we extract 
from Eq. |l9| for the behavior of the superfluid density 
with x are fairly independent of the precise functional 
forms of Wq and Ao. The only input that is needed is 
the existence of DDW order, with diminishing strength 
with x, and complementary development of the DSC or- 
der. The DDW order eats away part of the superfluid 
density from an otherwise pure DSC system even within 
the superconducting dome in the T — x phase diagram. 
With increasing x, the DDW order weakens, hence the 
superfluid density increases in a model-independent man- 
ner. 

In Fig. ^ we plot n s /m* against the doping-fraction x 



using Eq. |l5| and Eq. [TJ for Wo and A . To illustrate 
the qualitative robustness of our result, we also show, 
by the dashed curve, the same plot for Wo dispersing 
quadratically with x between the same two end points 
as in Eq. As can be seen from the figure the zero- 
temperature superfluid density, in a model-independent 
manner, shows a rapid drop in the underdoped regime, 
similar to that observed in experiments. 




FIG. 3. The T = superfluid stiffness n 3 /m* as a function 
of doping fraction x. The vertical axis is normalized by the 
superfluid stiffness at x — 0.2. Here, we have set t = 0.5 eV 
and the chemical potential is taken to be fi — —x 2 . The am- 
plitudes of the order parameters are as given in Eqs. [TB] and 
|l~6| . The dashed curve is n s /m* for Wo dispersing quadrati- 
cally with x as explained in the text. 

The leading temperature dependence of the superfluid 
stiffness can be evaluated from Eq. At temperatures 
much smaller than the relevant energy scales Wo and Ao, 
only the nodal regions close to the points (tt/2, tt/2) and 
symmetry-related points will contribute in the suppres- 
sion of the superfluid weight. By expanding around those 
points, one can see that the leading temperature depen- 
dence is indeed linear for the optimally and moderately 
doped samples in a fairly wide range of temperatures. 
For these doping concentrations, where A is larger than 
or comparable to W , W plays a subleading role to A 
in determining the temperature dependence of the sup- 
pression. This is due to the peculiar band structure of 
the problem. 

On the other hand, for the heavily underdoped sam- 
ples the situation is quite different. Critical fluctuations 
are clearly important to determine any low-temperature 
property of the system, for it is close to a quantum crit- 
ical point. Even if we ignore these fluctuations, and 
rely strictly on our mean-field results, the conclusions are 
very different from moderate or optimally doped samples. 
Though in the asymptotically low temperature regime, 
the depletion of the superfluid density is linear in tem- 
perature, there is an intermediate temperature range over 
which the suppression of the superfluid density actually 
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behaves as VT. As the DDW gap is much larger than the 
superconducting gap in these heavily underdoped sam- 
ples, Wo crosses over to produce the leading contributions 
in the expansion of Eq. around the nodes, and is even- 
tually responsible for the VT behavior of the suppres- 
sion. All of this is summarized in Fig. ||, where we have 
plotted the temperature-dependent part of the superfluid 
stiffness with temperature scaled by the superconducting 
gap, for six values of the doping fraction x. The plot 
shows VT behavior for the heavily underdoped systems 
for the experimentally relevant temperatures For 
asymptotically low temperatures for these doping values, 
and for a fairly wide range of temperatures for moder- 
ate or optimal doping, the data show exact collapse on 
a single straight line signifying a unique low-temperaure 
slope. 
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FIG. 4. The temperature-dependent part of the superfluid 
stiffness plotted against scaled temperature for six values of 
the doping. The doping concentrations are indicated in the 
legends. The values of Wo and Ao for a given doping are 
derived from Eqs. [IBJ and [l6|. /i is derived from /i = — x 2 and 
t = 0.5 eV. The dotted li ne is an extrapolation of the results 
to zero temperature. 



V. DISCUSSION 

The cuprate phase diagram proposed in reference 
is reproduced here in Fig. [l] for easy reference; we only 
show the relevant portion of the phase diagram, ignoring 
the complex set of competing charge ordered states in 
the underdoped regime as well as the spin glass phase. 

From our calculation, the behavior of the T — su- 
perfluid density with x is as expected from experiments 
in the underdoped regime of the high-T c cuprates. The 
striking rapid suppression of the superfluid density (and 
sharp increase of the London penetration depth) below 
x ~ 0.2 is naturally explained by the emergence of the 
singlet DDW at that doping-fraction. The behavior of 



the superfluid density near T c or near x « 0.05 (the low- 
est doping at which superconductivity occurs) will be 
strongly affected by thermal or quantum phase fluctua- 
tions, so our Hartree-Fock results will be suspect in those 
regimes. However, away from these critical points, we 
might expect our calculation to be on solid footing, and 
it is encouraging to see that it agrees with experimental 
measurements of n s . The temperature dependence also 
captures the striking linearity shown in Fig. |^, and the 
doping independence of the slope, along with a y/T de- 
pendence in an intermediate temperature regime in the 
heavily underdoped regime akin to the experiments [|l2| . 

The neutron scattering peaks observed in the pseudo- 
gap and superconducting regimes are broadly consistent 
with our calculations. There are peaks at (7r,7r) at ener- 
gies which are controlled by the single-particle gap and 
the doping. There are no observable peaks at (tt, 0) and 
symmetry-related points because, from our analysis, they 
would be at high energies controlled by the band disper- 
sion where they are likely to be strongly damped. This 
resolves an earlier puzzle of the analysis of the staggered 
flux phase Though we have not presented here re- 
sults for finite tempeartures, from the phase diagram and 
our zero temperature theory, the qualitative aspects at 
finite temperatures can still be explained. However, the 
details of the peak position in energy and its doping de- 
pendence is beyond a simple Hartree-Fock calculation. 

If one approaches along the path labeled (1) in Fig. [I], 
at first one will find no structure in the spin-spin corre- 
lator, typical of a Fermi liquid. Below T* , the correlator 
peaks at an energy 2Wq, twice the maximum of the DDW 
gap at the wavevector Q. This is consistent with experi- 
ments. Below T c , the intensity in the spin-spin correlator 
is amplified as the DSC order develops, as discussed in 
Ref. || . Due to the coexistence of DDW and DSC, the 
peak will be at a higher energy — 2\jW% + Aq, shifted 
by an amount that depends on the doping. To see this, re- 
call that, for coexisting DDW and DSC, the peak-energy 
shifts to higher values with which increases with dop- 
ing. 

Consider now the path (2) in Fig. [[]. In the DDW 
state, and for small chemical potential, the neutron scat- 
tering intensity should exhibit a peak at wavevector Q 
and energy 2Wq; recall that, for DDW order, the peak- 
energy shows no movement with \fi\, but is destroyed by 
a large \fi\ > Wo. Since the DDW gap varies with x 
as in Eq. [Ij, the peak would shift to higher energies as 
x is decreased, if xo alone were responsible for the spin- 
fluctuation spectrum. However, the reduced Coulomb in- 
teraction U is important since it shifts the peaks to lower 
energies, broadening them simultaneously. For smaller 
x, due to the lower density of the holes, the effective U 
increases, and, as a result, as shown in Fig. ^ the peak 
moves to lower energies. Consequently, the peak-energies 
in this part of the phase diagram would be influenced by 
these two competing effects. 
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Consider the path (3) in Fig. [[| In the coexisting re- 
gion, DDW plus DSC, the peak energy is 2y/Wg + A|, 
shifted by \fj,\, but now, in contrast to T > T c , the mag- 
nitude of the total order parameter remains constant as 
per Eq. [l^. Of course, U must be important as well, 
modifying this conclusion. 



ACKNOWLEDGMENTS 

We would like to thank Qiang-Hua Wang for pointing 
out a minor error in an earlier version of this paper. S. 
C. and S. T. are supported by NSF-DMR-9971138, and 
C. N. by NSF-DMR-9983544 and the A.P. Sloan Foun- 
dation. The work of H. -Y. K. was conducted under the 
auspices of the DOE, supported by funds provided by the 
University of California for the conduct of discretionary 
research by Los Alamos National Laboratory. 



[1] 
[2] 
[3] 
M 



[5] 



[6] 
[7] 



S. Chakravarty, R. B. Laughlin, D. K. Morr, and C. 

Nayak, Phys. Rev. B 63, 94503 (2001). 

C. Nayak, Phys. Rev. B 62, 488 (2000); ibid, R61 35 



(2000); C. Nayak and F. Wilczek, |cond-mat/9510132 
S. Chakravarty and H.-Y. Kee, Phys. Rev. B 61, 
14821(2000). 



H. A. Mook, P. Dai, F. Dogan, |cond-mat/0102047| ; H. A. 
Mook, Pengcheng Dai, S. M. Hayden, A. Hiess, and F. 
Dogan, "Polarized neutron beam observation of antifer- 
romagnetic order in YBa2Cu306.e" , preprint. The data 
of Y. Sidis et al., Phys. Rev. Lett. 86, 4100 (2001), on 
lower doping material YBCO6.5 are entirely different and 
appear to reflect remnant spin magnetism instead. There 
is no contradiction in this, because DDW order can dis- 
appear as the doping is reduced; see S. Ch akravarty, H. 
-Y. Kee and C. Nayak, |:ond-mat/0101204 . 
J. E. Sonier, J. H. Brewer, R. F. Kiefl, R. I. Miller, G. D. 
Morris, C. E. Stronach, J. S. Garner, S. R. Dunsiger, D. 
A. Bonn, W. N. Hardy, R. Liang, R. H. Heffner, Science, 
to be published. 



M. V. Eremin et al, |cond-mat/0105508 . 

For recent reviews, sec H. F. Fong et al, Phys. Rev. B 

61, 14773 (2000); P. Dai et al, Phys. Rev. B 63, 054525 

(2001). 

For a cross section of theoretical papers on the neutron 
resonance peak, see K. Maki and H. Won, Phys. Rev. 
Lett. 72, 1758 (1994); M. J. Lercher and J. M. Wheatley, 
Phys. Rev. Lett. B 49, 736 (1994); M. Lavagna ang G. 
Stenmann, Phys. Rev. B 49, 4235 (1994); G. Stenmann 
and C. Pepin, and M. Lavagna, Phys. Rev B 50, 4075 
(1994); F. Onufrievaand J. Rossat-Mignod, Phys. Rev. B 
52, 7572 (1995); G. Blumberg, B. P. Stojkovic, and M. V. 
Klein, Phys. Rev. B 52, 15 741 (1992); A. A. Abrikosov, 
Phys. Rev. B 57, 8656 (1998); L. Yin, S. Chakravarty, 
and P. W. Anderson, Phys. Rev. Lett. 78, 3559 (1997); 



N. Bulut and D. J. Scalapino, Phys. Rev. B 53, 5149 
(1996); D. Z. Liu, Y. Zha, and K. Levin, Phys. Rev. 
Lett. 75, 4130 (1995); I. I. Mazin and V. M. Yakovenko, 
Phys. Rev. Lett. 75, 4134 (1995); J. Brinckmann and 
P. A. Lee, Phys. Rev. Lett. 82, 2915 (1999); Y. J. Kao, 
Q. Si, and K. Levin, Phys. B 61, 11 898 (2000); D. K. 
Morr and D. Pines, Phys. Rev. Lett. 81, 1086 (1998); 
A. J. Millis and H. Monien, Phys. Rev. B 54, 16 172 
(1996); T. Tanamoto, H. Kohno, and H. Fukuyama, J. 
Phys. Soc. Jpn. 61, 1886 (1992); B. Normand, H. Kohno, 
and H. Fukuyama, J. Phys. Soc. Jpn. 64, 3903 (1995); 
A. Abanov and A. V. Chubukov, Phys. Rev. Lett. 83, 
1652 (1999); S. Sachdev, C. Buragohain, and M. Vojta, 
Science 286, 2479 (1999); E. Dernier and S. C. Zhang, 
Phys. Rev. Lett. 75, 4126 (1995); E. Demler, H. Kohno, 
and S. C. Zhang, Phys. Rev. B 58, 5719 (1998). 
[9] Pengcheng Dai, H. A. Mook, G. Aeppli, S. M. Hayden, 
and F. Dogan, Nature 406, 965 (2000). 

[10] JL. Tallon and J.W . Loram, Physica C 349, 53 (2001), 
|cond-mat/0005463| ; J.W.Loram et al, J. Phys. Chem. 
Solids 59,2091 (1998); J.L.Tallon et al, Phys. Stat. Sol. 
B 215, 531 (1999). 

[11] Y. J. Uemura et al, Phys. Rev. Lett, 62, 2317 (1989). 

[12] C. Panagopoulos, J. L. Tallon, and T. Xiang, Phys. Rev. 
B 5 9, R6635 (1999); see also D. Brown, Ph. D. thesis, 
Cambridge University, unpublished (2000). 

[13] E. Cappelluti and R. Zeyher, Phys. Rev. B 59, 6475 
(1999); E. Cappelluti and R. Zeyher, Physica C 312, 313 
(1999); E. Cappelluti and R. Zeyher, Europhys. Lett. 49, 
487 (2000); R. Zeyher, Physica C 341,121 (2000). 

[14] I. Affleck and J. B. Marston , Phys. Rev. B, 37, 3774 
(1988); J. B. Marston and I. Affleck, Phys. Rev. B, 39, 
11538 (1989). 

[15] E. Fradkin, Field theories of condensed matter systems 

(Addison- Wesley, Redwood City,1991). 
[16] A. Ino et al, Phys. Rev. Lett, 79, 2101 (1997) 
[17] N. Furukawa and M. Imada, J. Phys. Soc. Jpn. 62, 2557 

(1993) 

[18] E. Dagotto et al, Phys. Rev. Lett, 67, 1918 (1991). 

[19] See, for example, S. A. Kivelson and V. J. Emery, 
in "Strongly Correlated Electronic Materials: The Los 
Alamos Syposium 1993," edited by K. S. Bedell et al 
(Addison Wesley, Reading, 1994) p. 619. 

[20] V. J. Emery, S. A. Kivelson, and O. Zachar, Phys. Rev. 
B 56, 6120, (1997). 

[21] J.R.Schrieffer, The Theory of Superconductivity (W. 
A. Benjamin, Inc., Reading, 1964). We fold the usual 
Nambu-basis, which is defined in the full Brillouin zone, 
onto the reduced Brillouin zone to incorporate DDW si- 
multaneously with DSC; see also F. C. Zhang, Phys. Rev. 
Lett. 64, 974 (1990) 

[22] See also a re lated preprint, Q. -H. Wang et al, cond 



mat /00 11398 



[23] J. W. Loram et al, Physica C, Feb 20-25, 2000. 
[24] D. J. Scalapino et al, Phys. Rev. B, 47, 7995, (1993). 
[25] P. Miller et al, Physica C, 210, 343 (1993); G. Bilbro 
and W. L. McMillan, Phys. Rev. B, 14, 1887 (1976). 



8 



